Shapes of hydrophobic thick membranes 
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We introduce and study the behavior of a tethered membrane of non-zero thickness embedded 
in three dimensions subject to an effective self-attraction induced by hydrophobicity arising from 
the tendency to minimize the area exposed to a solvent. The phase behavior and the nature of 
the folded conformations are found to be quite distinct in the small and large solvent size regimes. 
We demonstrate spontaneous symmetry-breaking with the membrane folding along a preferential 
axis, when the solvent molecules are small compared to the membrane thickness. For large solvent 
molecule size, a local crinkling mechanism effectively shields the membrane from the solvent, even 
in relatively flat conformations. We discuss the binding/unbinding transition of a membrane to a 
wall that serves to shield the membrane from the solvent. 



INTRODUCTION 

The structure and dynamics of sheets or surfaces are 
important in diverse contexts in everyday life, in cell bi- 
ology, in field theory and quantum gravity, and in con- 
densed matter physics 0, Q . Previous studies have con- 
centrated on idealized surfaces - the thickness is neglected 
principally to facilitate an analytical approach (see the 
excellent work by David, Duplantier, Guitter and Wiese, 
reviewed in Ref . [3] ) or partially captured by tethering to- 
gether, in a two dimensional array, hard spheres, whose 
diameter equals the surface thickness. In the latter case, 
both the steric and interactions induced by the environ- 
ment or solvent are modeled as standard two-body po- 
tentials which become singular (Dirac delta function) in- 
teractions in the continuum limit Q. In this limit, self- 
interacting manifolds, in 1 dimension corresponding to a 
polymer chain and in 2 dimensions corresponding to a 
membrane, the subject of our study, cannot be correctly 
described by two body interactions The difficulty 
arises because a pair-wise interaction can only depend on 
the distance between the two parts of the manifold. One 
is unable to account for whether these two parts are near 
each other just because they lie on adjacent parts of the 
manifold or whether they are distant along the manifold 
and thence truly interacting. One could cure this prob- 
lem by introducing some microscopic cut-off such that 
nearby parts of the manifold, at a distance smaller than 
this cut-off along the manifold, are defined to be non- 
interacting. This procedure is expected not to affect the 
long distance behavior but it affects many important in- 
termediate scale effects of interest. 

Extensive studies 0-0| have shown that there are 
vastly different behaviors between a conventional poly- 
mer modeled as a chain of non-overlapping spheres and 
a tube due distinct symmetries in the two cases. The 
latter can be imagined as a chain of coins whose symme- 
try is cylindrical. Self-avoidance of a flexible tube can 
be imposed by using a three-body potential A tube 



subject to maximal compaction takes on the geometry 
of an optimal helix [Bj], whose local curvature is equal to 
the tube thickness. Strikingly, this geometry is adopted 
by a-helices in proteins. On the other hand, a fee lattice 
is an optimal arrangement for a chain of spheres sub- 
ject to compaction, provided the tether constraints do 
not conflict with this arrangement. Surprisingly, chains 
of overlapping spheres Q behave in a manner akin to 
a tube. This is because the overlap between adjacent 
spheres along the chain result in a loss of spherical sym- 
metry and endow an anisotropy to the chain much as a 
tube does [Icj ]. 

The principal theme of this paper is to introduce a 
model of a thick hydrophobic membrane that encapsu- 
lates the correct symmetry associated with its intrinsic 
two dimensional nature. The implementation of the non- 
zero thickness of the membrane is carried out through 
a four-body potential Unlike the model of tethered 
spheres, the membrane thickness in our model provides 
a length scale that is independent of the discretization 
of the membrane surface, and thus is free of the sin- 
gularity problem in the continuum limit. The inter- 
action with the solvent molecules is captured through 



Hadwiger's theorem [ll| as explained in [121] : the ef- 
fective self-interactions of the membrane, induced by a 
short-range solvent-membrane potential, can have only 
four distinct contributions. Two of these interactions are 
proportional to the integrals of the mean curvature and 
Gaussian curvature over the entire surface. The other 
two interactions pertain to an entropic contribution aris- 
ing from the volume excluded to the solvent molecules 
[13| and an energy contribution measured by the area of 
the membrane that is exposed to the solvent molecules. 
The membrane thickness already takes into account the 
mean curvature effects because, at each position on the 
surface, the curvature radii have to be larger than the 
thickness itself. The integral of the Gaussian curvature 
is constant if the membrane topology is kept fixed. Fur- 
thermore, because the membranes we consider are open, 
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the volume excluded to the solvent molecules is partly 
taken into account by the buried area or equivalently by 
its complement, the exposed area. Thus our focus here is 
on understanding how the area of the membrane exposed 
to the solvent molecules influences the nature of the con- 
formations of a hydrophobic membrane immersed in a 
solvent. We show that distinct behaviors are obtained 
for different ratios between the membrane thickness and 
the solvent molecule size, both in the absence and pres- 
ence of an inert wall, which serves to shield one side of 
the membrane from the solvent. We find several familiar 
conformations that have not been observed in previous 
studies of conventional models. 



MATHEMATICAL DESCRIPTION OF THICK 
MEMBRANE 

The notion of thickness can be introduced for a contin- 
uum surface in 3-D. The position vector of a point on the 
surface is given by r(s), where s— (s 1 ,,? 2 ) G D defines 
the local (curvilinear) coordinate frame on the surface. 
The two tangent vectors are 



(1) 



The square of the distance between two neighboring po- 
sitions on the surface, s, s + ds is given by 



where 



dr 2 = gjjdsMs 3 , 



gij = dir(s) ■ djr(s) 



(2) 



(3) 



is the metric tensor. We will assume that the tangent 
vectors are linearly independent at all points on the sur- 
face which is equivalent to 



[d 1 r(s)xd 2 T(s)j = det( 9ij ) = g ? . (4) 

The surface element at position s is given by 

dS(s) = d 2 s ^g~[g) . (5) 

We also assume that the local coordinate frame is chosen 
such that the normal to the surface 



h(s) = 



dir(s) x d 2 r(s) 
\dir(s) x d 2 r(s)\ 



(6) 



is always on one side of the surface. If the thickness of the 
surface is 2A, then the radius of the sphere tangent at 
r(s) and going through another point r(t) on the surface 
(t is a vector in the curvilinear coordinates), 



K = 



f(t)-r(s)] 2 



2n( a >[r(t)-r(§)] ' 



(7) 



cannot be smaller than A for all s and t (fig. 1). 

For a triangulated surface with fixed connectivity cor- 
responding to the middle surface of a membrane, assume 
that the vertices are denoted as i, j, k . . . whose positions 

are given by rj, r^, A triangle of vertices i, j, and 

k is denoted by (ijk). For each triangle (ijk), we assign 
an orientation such that for the completely fiat surface 
the normal vector 



(r, - rj) x (rfc -rj) 

l( r J - r X ( r k -Tj)\ 



(8) 



always points in the same, say up, direction. The thick- 
ness is defined in analogy with the continuum case. The 
radius of the sphere tangent to the triangle (ijk) at its the 
center of mass r^ = (r$ + rj + rfc)/3 and going through 
the node r m , 
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ijkm — 



(9) 



is constrained to be greater than or equal to 
A V(ijk),m^ i,j,k. 



EXPOSED AREA OF THICK MEMBRANE 

Assume that the solvent molecule is a spherical paint 
brush of radius R (fig. 1). The exposed area of a mem- 
brane is defined as the membrane's surface area accessible 
by this spherical paint brush. Each membrane has two 
sides, + and — , defined by two triangulated surfaces, S + 
and <S_, at a distance A from the middle surface. We 
assign a normal to each node r 4 of the middle surface as 
the average normal of all the triangles sharing the z-th 
node: 



(10) 



where the sum is over all oriented triangles sharing the 
i-th node. The number of such triangles is six in the 
bulk and less than six at the boundary. Each n i; by 
definition, is oriented in the direction of S + . Thus we 
define the nodes of the surfaces S± as 



r, : ± Am. 



(11) 



To calculate the exposed area of the surfaces S±, one 
places a sphere of radius R tangentially to each triangle 
on the surfaces S± at its center of mass, and then to 
check whether the sphere overlaps with any of the other 
points on these two surfaces. If there is no overlap the 
triangle is said to be exposed to the solvent. The follow- 
ing quantity can be calculated for each triangle (ijk) in 
the <S± surfaces: 
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(12) 
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FIG. 1: Cartoon of a thick membrane and a solvent molecule. 
Each side of the membrane is at a distance A from the mid- 
dle surface (dashed line). Self-avoidance requires that any 
sphere tangent to the middle surface at one point and pass- 
ing through another point on this surface must have a radius 
r > A. The solvent molecule is considered as a spherical 
paint brush of radius R. The membrane exposed surface area 
is defined as the area on both side of the membrane that is 
accessible to that spherical paint brush. 



where p,^ fe and N^- fe is the center of mass and the normal 
of that triangle, respectively, and A can be + or — . The 
exposed areas in the surfaces S± are defined as 
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(13) 



where 0() is the step function, and the sum is taken over 
all triangles on a given surface. 



MODEL OF THICK HYDROPHOBIC 
MEMBRANE 

Following standard practice, here we model the thick 
membrane as a set of points forming a triangular lattice 
representing the membrane's middle surface. The links 
on the lattice are fixed and have lengths that are allowed 
to vary freely between 0.9/ and 1.1/, where / is the unit 
length. Self-avoidance is implemented through an effec- 
tive four-body potential (Eq.(9)) to ensure the membrane 
thickness of 2A. Due to the effects of discretization, this 
procedure is accurate as long as A// > 0.5. Note that 
A provides a length scale that is independent of the dis- 
cretization length /. The continuum limit can be ap- 
proached by decreasing / while keeping A constant. 

We consider a membrane that is hexagonal in shape 
when it is perfectly flat with the edge length of L. The 
solvent molecule is modeled as a spherical paint brush 
with radius R. The interaction between solvent molecules 
is not considered in our study and yields a higher order 
correction. For a hydrophobic solvent, there is a ten- 
dency to minimize the exposed area by appropriate fold- 
ing of the membrane, which, at non-zero temperatures, is 
counteracted by an entropic contribution promoting rel- 
atively flat conformations. An effective self- attraction is 



thus introduced through a potential energy: 

£ = (£/z 2 )(£ + + ir), 



(14) 



where e provides an energy unit and S are the exposed 
areas of the + and — sides, respectively. 

Monte Carlo simulations are carried out to study the 
behavior of the membrane at various temperatures. In 
a Monte Carlo move, one node in the triangulated lat- 
tice of the membrane's middle surface is selected and 
displaced in a random direction by a random displace- 
ment. The magnitude of the displacement is constrained 
to be less than 10% of the lattice constant /. The move 
is rejected immediately if (1) a new bond length on the 
lattice due the displacement is smaller than 0.9/ or larger 
than 1.1/, or (2) the membrane thickness constraint is vi- 
olated. Otherwise, the move is accepted with probability 
P = min[l, exp(— AE/ksT)], where T is temperature, 
and AE is a change in the effective energy. A parallel 
tempering scheme with 16 to 20 replicas is adopted to 
get efficient sampling of the conformational space. 



LARGE VS. SMALL SOLVENT SIZE 

Figs. [2] show typical conformations of a membrane in 
the large and small solvent size regimes. At high tem- 
peratures, the entropy dominates, there is no tendency 
for minimizing the area exposed to the solvent, and the 
conformations in both regimes are flat. On lowering the 
temperature, the membranes tend to become compact by 
minimizing their exposed area to the solvent but their 
behaviors are remarkably different for the two cases of 
small and large solvent molecule sizes. The correlations 
between the normal vectors of the membranes also show 
a major difference between the low temperature confor- 
mations in the large and small solvent regimes (Fig. [3|). 

When R is larger than A, the low temperature phase 
is crinkled and contains many different conformations. 
The crinkling is essentially a local deformation that does 
not violate the local bending constraint of a thick mem- 
brane yet excludes the large paint brush from accessing 
the surface. Locally, below the persistence length, the 
membrane is essentially flat yet crinkled. At the low- 
est temperature, the crinkled conformations are almost 
completely buried (Fig. [2^,). We have verified that, for 
the same solvent-induced energy function, the impenetra- 
ble plaquette model [3], behaves similarly to the thick 
membrane in the large solvent molecule limit - i.e. the 
thickness is not relevant when the solvent molecules are 
large. Interestingly, for large solvent size, the local crin- 
kling mechanism induces a stiffening of the membrane in 
the flat phase (Fig. [3]). 

When R is considerably smaller than A, the mem- 
brane exhibits a transition from the high temperature flat 
phase to the low temperature folded phase characterized 
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(a) (b) 




R=0.5A , T=0.2 R=0.5A , T=0.1 



FIG. 2: Conformations of a hydrophobic membrane im- 
mersed in water. The membrane has an edge length of L — 71 
and the thickness 2A = 1.6 1, (a) Two views of a conforma- 
tion found at T = O.le/fcs for a membrane when the solvent 
radius, R — 5A, is large. The triangular lattice shown indi- 
cates the membrane's middle layer. The membrane is fully 
buried due to the local crinkling mechanism. Note the the 
membrane surface is rough but the overall shape is flat, (b- 
d) Three typical folded conformations at low temperatures for 
a membrane in the small solvent regime. The solvent radius is 
R = A/2, and the temperatures are T — QAe/ks (b), 0.2e/fc_e 
(c), and O.le/fcs (d) as indicated. The buried fractions are 
0.31, 0.37 and 0.44 for the three cases (b, c, d), respectively. 
Note that folding occurs predominantly along one axis. 
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FIG. 3: Correlation between normal vectors. Depen- 
dence of the correlation between normal vectors, {n^riaj+r) — 
{n x ){n x + r ), on r, the distance between them, are shown for 
the large solvent (7? = 5A) and small solvent (R = 0.5A) 
cases. The small solvent case is shown as inset figure. The 
membrane has the same characteristics as the one shown in 
Fig. [2] The correlation functions are shown for two different 
temperatures: very high, T — 100e/fcs (solid line) and low, 
T = O.le/fcB (dotted line), as indicated. The data shown for 
each temperature are calculated from 100 statistically inde- 
pendent conformations sampled in parallel tempering Monte 
Carlo simulations under equilibrium conditions. The oscil- 
lation in the correlation function for small solvent size at 
low temperature reflects the folding along a preferential axis. 
Note that in the large solvent case, the correlation, at fixed 
distance, is larger as temperature decreases demonstrating a 
stiffening of the membrane due to local crinkling. 



by a few different folded geometries (Fig. [3b,c,d). The 
membrane undergoes a spontaneous symmetry breaking 
within the plane of the membrane. The folding of the 
membrane occurs multiple times along a dominant fold- 
ing axis. 

RIBBONS IN 2-D 

We have carried out also simulations of solvophobic 
ribbons in 2-D (Fig. Q|. The ribbons have a thickness 
of 2A and the solvent molecule in this case is considered 
as circular disc of radius R. We consider a discretized 
ribbon axis described by a set of equally spaced points 
separated by b = A/2. The self-avoidance constraint is 
imposed by requiring that none of the radii of the circles 
passing through each of the triplets of the points on the 
ribbon axis is smaller than A. All triplets of points are 
considered. Like for membranes, we found two regimes 
related to the solvent size. For R > A, the ground state 
of the ribbon is a straight conformation with local crin- 
kling (Fig. @Ji). For R < A, one finds folded conforma- 
tions with some local crinkling on the outer perimeter 
(Fig. 2] b,c). The folded conformations of the ribbon 



correspond to a folded membrane in a dimcnsionally re- 
duced representation. For sufficiently small solvent size, 
the ground state of the system is akin to a rolled carpet 
(Fig. St). 

Using simple analytic calculations, we find that when 
R > A(2— 6 2 /A 2 ), a ribbon can adopt a virtually straight 
conformation with local crinkling that buries the perime- 
ter completely. Such a conformation is shown in Fig. 0^, 
for R = 2 A (the buried fraction does not reach 1 only 
because of edge effects). In an optimal local crinkled 
conformation, the radius of curvature is equal to A at 
all points. In the continuum, one can approximate a 
crinkled conformation of the ribbon axis as being formed 
by joining semicircles of radius A, the first convex up- 
wards, the second concave upwards, and so on. The 
buried fraction of such a continuum crinkled ribbon is 
given by / = 1 — {2 /it) arcsin ( 2 +_r/a ) ' wm ch converges 
to 1 as R/A — > oo. 





R=A 



R=0.25A 



FIG. 4: Conformations of a solvophobic ribbon in 2-D in the 
presence of a solvent. In the regime where the membrane is 
folded (Fig. [2} this is a simpler (dimensionally reduced) repre- 
sentation. The ribbons shown are of thickness 2A and length 
L — 25 A (a) and L = 75A (b, c). The solvent is comprised 
of circular discs of radius R whose value is indicated. The 
ribbon axes are shown as dashed lines. The ribbon perime- 
ters are shown in thinner (thicker) solid line depending on 
whether they are buried from (exposed to) the solvent. The 
buried fractions are 0.94 (a), 0.82 (b), and 0.72 (c) for the 
three cases shown. Conformations (b) and (c) are the ground 
states obtained in Monte Carlo simulations that seek to min- 
imize the exposed perimeter of the ribbons. Conformation 
(c) corresponds to the familiar rolled carpet conformation of 
membranes in 3-D in which the normal vectors lie isotropi- 
cally within the plane perpendicular to the folding axis. 



MEMBRANE NEAR A WALL 



We have also studied the behavior of a tethered mem- 
brane of non-zero thickness near a wall which provides 
a new mechanism for burying the area of a membrane 
adsorbed on it. For small solvent molecule size, there 
is a competition between the binding to the wall and 
folding for reducing the exposed area. Our simulations 
show that below a temperature, which increases with 
membrane size, and is higher than the folding temper- 
ature, the membrane is adsorbed on the wall (Fig. [5]). 
An adsorbed membrane would be flat at intermediate 
temperatures and folded at lower temperatures. Similar 
considerations suggest that multiple hydrophobic mem- 
branes would adhere to each other, at any temperature, 
to shield much of the surface from the solvent yielding 
an effectively rigid assemblage of membranes. For large 
solvent molecules, the adsorption happens below a tem- 
perature, that is higher than the temperature of effective 
local crinkling. At very low temperatures, the membrane 
is crinkled while remaining bound to the wall. 



T=3.5 
T=2.5 
T=1.5 



(b) 30 

20 
10 

& 

Ll_ 

-10 
-20 
-30 



FIG. 5: Adsorption of a membrane on the wall, (a) Tem- 
perature dependence of the specific heat for the system of a 
hydrophobic membrane and a wall. The membrane is hexag- 
onal in shape and has the edge length L = 51 and the thick- 
ness 2 A = 1.2L The solvent radius is R — A/2. Adsorp- 
tion is induced by the shielding of membrane surface from 
the solvent when it is found near the wall. The adsorption 
transition is observed at temperature T « 2.5e/fcs. Mem- 
brane conformations are shown for selected temperatures as 
indicated. Both the bound and unbound states at interme- 
diate and high temperatures are flat. The ground state how- 
ever is a folded membrane bound the the wall, and can be 
found only at very low temperatures (T = O.le/fcs). (b) 
The dependence of the effective free energy, F, on the dis- 
tance, z cm , of the membrane center of mass to the wall for 
three temperatures: above, below and at the adsorption tran- 
sition as indicated. The effective free energy is calculated 
as F(T,z cm ) = -k B TlnP(T,z crn ), where P(T,z cm ) is the 
probability of finding the membrane at the distance z cm of its 
center of mass to the wall at temperature T. 



DISCUSSION 

Previous computational studies of membranes have 
considered, in essence, models in the infinitesimal thick- 
ness limit or 2-d array of spherical particles. The key 
finding of studies of traditional models (for a review see 
Ref.[3|) is that a self- avoiding tethered membrane always 
exists in a flat phase even without any bending rigid- 
ity [I34la |. The inclusion of self- attraction 17 1 leads to 
several folded phases at intermediate temperatures along 
with hints of a crumpled phase 18[ and an isotropic 
collapsed phase at low temperatures. A recent study 
[l9l | with a Lennard- Jones potential yielded folded sheets 
and cylindrical conformations. The crumpled state has 
been also studied for an elastic sheet confined in an im- 
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penetrable sphere [§c| or folded under external forces 
21]. Studies of self- intersecting membranes (or phantom 
membranes) predict a crumpled phase in which the mem- 
brane's radius of gyration, R g , grows as y/lnL, where 
L is the membrane's linear size 22J, |23[. With an in- 
crease in bending rigidity, phantom membranes undergo 
a second-order transition from the crumpled phase to 
the flat phase [3, l2S| . The flat, crumpled, and col- 
lapsed phases have been observed experimentally in red 
cell membrane skeletons [26j j . graphite oxide membrane 
27 1, and molybdenum disulfide [281 ] . 

Here we show that incorporating non-zero thickness 
yields qualitatively new behavior of a hydrophobic mem- 
brane seeking to minimize the area exposed to the sur- 
rounding solvent molecules. Our simulations show that 
in the small solvent molecule regime, the thickness of the 
membrane introduces a spontaneous symmetry breaking 
within the plane of the membrane. The folding of the 
membrane occurs multiple times along a dominant fold- 
ing axis. This is akin to the everyday experience of the 
greater ease of folding a sheet of paper in a fluted man- 
ner compared to first folding the paper along one axis 
followed by folding the now two layered sheet along a per- 
pendicular axis and so on. A rolled carpet conformation 
is also effective in burying much of the surface from the 
solvent molecules (Fig. 2b). Note that the folded struc- 
tures obtained here are similar to the those in the tubule 
phase [29I ] - even in the absence of any anisotropic bend- 
ing energy, we observe spontaneous symmetry breaking. 
For a large membrane, one obtains uncoordinated folding 
at the edges into metastable conformations. The struc- 
tures of the folded thick membrane are qualitatively sim- 
ilar to those observed in molybdenum disulfide films [28[ 
and in folding of viscous sheets and filaments |Ioj ■ 

For large solvent molecules, the local crinkling mech- 
anism effectively shields the surfaces of the membrane 
from the solvent, while preserving the overall flat topol- 
ogy. Similar local deformations are observed in real mem- 
branes. Wrinkled rigid structures have been observed 
experimentally upon cooling partially polymerized phos- 
pholipid membranes 31( and have been interpreted as 
resulting from quenched curvature disorder leading to a 
glassy phase. Studies of wrinkled patterns of polymer 
films [3 21 ] have indicated that thickness plays a role in 
the formation of these patterns. 

We note that the our study considers paint brush sizes 
in two distinct regimes: smaller than or larger than the 
membrane thickness. Membranes thicker than the size 
of a water molecule (« 2. 8 A) are quite common. The 
other limit is obtained when the role of solvent is played 
by large complexes such as proteins (3-10 ran), t-RNA (7 
nm), antibodies (12 nm) and ribosomes (20-30 nm). In- 
terestingly, for a relatively dilute concentration of these 
complexes, even of the order of 10% in volume fraction 
according to Asakura-Oosawa theory (33[, the entropy 
depletion effect alone is sufficient to create a tendency 



to fold the membrane leading to minimization of the ex- 
posed area. 

We are indebted to Mehran Kardar, Tom Luben- 
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manuscript. This article was supported by the Viet- 
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